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Abstract 

We argue that the measurement of open charm gives a unique opportunity to test the 
validity of pQCD-based and statistical models of nucleus-nucleus collisions at high 
energies. We show that various approaches used to estimate D-meson multiplicity 
in central Pb+Pb collisions at 158 A-GeV give predictions which differ by more 
than a factor of 100. Finally we demonstrate that decisive experimental results 
concerning the open charm yield in A+A collisions can be obtained using data of 
the NA49 experiment at the CERN SPS. 
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1 Introduction 



The measurement of production of D-mesons in nucleus-nucleus collisions at CERN SPS 
energies is a challenging experimental adventure. This is due to the short life time and 
expected low multiplicity of D-mesons which cause the small experimental signal origi- 
nating from D decays to be hidden in the large background composed of combinations of 
'non-signal' tracks. 

There is, however, increasing interest in results on charm production in nucleus- 
nucleus (A+A) collisions which motivates various experimental groups to design upgrades 
of the existing experiments or to build new experiments, which should allow to measure 
open charm production in A+A collisions 0, |^ . 

The aim of this paper is to demonstrate that the first measurement of D-meson 
production may be even possible using the current set-up of the NA49 experiment [Q. 

The paper starts with a summary of physics motivations for charm measurement in 
A+A collisions at the CERN SPS (Section 2). The expected multiplicity of D-mesons 
in central Pb+Pb collisions at 158 A-GeV is discussed in Section 3. An estimate of the 
statistical resolution of D-meson measurement which can be achieved using data of the 
NA49 experiment is given in Section 4. Summary and conclusions (Section 5) close the 
paper. 
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2 Physics Beyond QGP: pQCD vs Thermodynamics 



The main motivation for a broad experimental program in which nucleus-nucleus collisions 
at high energies are studied is the search for the Quark Gluon Plasma ||^, |^ . An impressive 
set of data was collected during the last decades and many unexpected phenomena were 
discovered 0]. However, the question whether a transient QGP state is created in the 
early stage of the collisions is still under vivid discussion. In our opinion this situation is 
due to the fact that there is no consensus on which models of high energy hadronic and 
nuclear collisions should be used to interpret experimental results. 

The data are in qualitative and quantitative agreement with the hypothesis of QGP 
creation in A+A collisions at the CERN SPS if a statistical model of the early stage 
is used for their interpretation [Q. This model is, however, considered a non-orthodox 
approach to hadronic and nuclear collisions as its basic assumptions can not be derived 
from the commonly accepted theory of strong interaction, QCD. Its validity is, therefore, 
more surprising than the conclusions reached within this model concerning QGP creation 
at the CERN SPS. 

On the other hand it is difficult to use QCD for the interpretation of the experimental 
results. Problems arise because almost all effects intuitively expected in the case of a 
transition to QGP are in the domain of so-called soft processes in which experimentally 
testable strict predictions of QCD are not possible. Attempts to build phenomenological 
QCD-inspired models of soft processes are not very succesful 0. Conclusive interpreta- 
tion of the data within these models seems to be impossible as one can not estimate the 
uncertainties introduced by the used approximations. 

As a possible solution one tries to identify phenomena sensitive to the early stage in 
which so-called hard QCD processes are believed to prevail |jlO|. In this case one hopes to 
obtain testable QCD predictions using perturbative methods. A well known example is 
the analysis of J /ip production in A+A collisions performed under the assumption that the 
initial production of cc pairs is described by pQCD. However the limits of the applicability 
of pQCD are theoretically not well defined. Therefore any assumption based on validity 
of pQCD calculations and used to interpret the data should be tested experimentally. In 
the case of the J / tp example this test can be done by the measurement of the open charm 
yield and its comparison to pQCD predictions. 

In connection with the above discussion of QCD-based models it is interesting to 
consider a hypothetical model in which, due to technical problems, strict experimentally 
testable predictions can not be obtained. The model gives approximate predictions, how- 
ever, the error made due to the approximations is difficult to estimate. Experimental 
results which are in agreement with the predictions of the model can be treated as a 
'proof of its validity. On the other hand any disagreement between the model and the 
data can be interpreted as due to the used approximations. Therefore the model can 
not be falsified |TT| before a substantial improvment in its predictive power. This ex- 
ample illustrates the logical problem which one may encounter discussing the validity of 
QCD-based models. 

In the summarized exciting situation concerning our understanding of strong interac- 
tions the role of the experimental results on open charm production in A+A collisions is 
unique. 



3 



Data on open charm should allow to test limits of applicability of the perturbative 
QCD methods. 

It is often assumed that the charm quark is heavy enough for pQCD treatment . 
Based on this hypothesis results on J/ ip production in A+A collisions are interpreted 
in the so-called 'suppression' models |T3|, |1^], in which the suppression is calculated 
relative to the J/tp yield expected from the pQCD factorization theorem and 
data on J/ip production in p+p interactions. It was, however, recently shown that 
the J ftp multiplicity is, in good approximation, proportional to pion multiplicity 
[^. Thus the J/ip yield shows an A-dependence which is characteristic of soft 
QCD processes. This experimental observation motivates the question whether 
pQCD can be used as a model of charm quark production in hadronic and nuclear 
collisions. 



• Data on open charm should allow to test the validity of a non-orthodox statistical 
model of the early stage of A+A collisions 

In this model charm quark multiplicity is large enough to justify the use of the 
thermodynamical approximation. Consequently not only the absolute yield of open 
charm but also its A-dependence are very different for the statistical and the pQCD- 
based models. 
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3 Multiplicity of D— Mesons in Nuclear Collisions 



This section starts with a review of the data on and D production in proton-nucleon 
interactions, which leads to an estimate of the mean multiphcity of + D mesons, 
(Z)° + D ), in nucleon-nucleon, N+N, interactions at 158 A-GeV. This result is further 
used to predict the multiplicity of D'^ + mesons in central Pb+Pb collisions at 158 
A-GeV using four different approaches. 



3.1 {D^ + D ) in N+N interactions 



There are six measurements of D-meson production cross section at various collision 
energies ranging from 200 GeV to 800 GeV, which allow to estimate the mean multiplicity 



of + in p+N interactions [T|, P, [T8|, |T|, |20|, HTJ. In p+p interactions a{D^ + D ) 
is measured at 360 GeV [T^, 400 GeV [jg and at 800 GeV At 200 GeV 0, 250 



GeV [|T7| and 800 GeV ||21| the cross section in p+N interactions is estimated based on 
measured cross section for p+A interactions assuming that (t(D° + Z)") is proportional 
to A. We note here that this assumption is not fully justified by the experimental results 
15 1 and may lead to an additional systematic error in estimating a{D^ + D ) in p+N 
interactions. In addition a symmetry of the xp distribution of D'^ + mesons with 
respect to xj? = is assumed for the calculation of the integrated cross section using 



measurements at Xi? > 0. In the case of the measurement at 200 GeV [0 only results 
for the sum of all D-mesons are published; in this case we assume that 50% of them are 
+ D mesons. 



^0 



The mean multiplicity of D + D mesons in p+p and p+N interactions is calculated 



as a ratio + D^)/a, where a is the inelastic cross section for p+p interaction at 

the corresponding energy. The values of {D^ + D are given in Table 1 and plotted 
as a function of \/s in Fig. 1. The results at ^/s = 19.4 GeV and y/s = 21.7 GeV 
differ by a factor of about 8, which suggests a possible large systematic uncertainty of 
the measurements at low energy not accounted for in the quoted errors. Taking this into 
account we estimate the mean multiplicity oi D° + D at ^ = 17.3 GeV {Elab = 158 
GeV) in nucleon-nucleon interactions by an arithmetic average of the measurements at 
^/s = 19.4 GeV and ^/s = 21.7 GeV, which gives: 



2.0- lo- 



ll 



A large systematic error up to several times the estimated value should be kept in mind. 



3.2 {D^ + D ) in Central Pb+Pb Collisions at 158 A GeV 

We discuss here four approaches which allow to estimate the mean multiplicity of + 
in central Pb+Pb collisions, {D^ + D )pbPb, at 158 A-GeV. It is assumed that the central 
collisions are selected by accepting only events with a low number of spectator nucleons 
from the projectile nucleus and that typical trigger conditions of the NA49 experiment 
are used [Q. This results in a mean number of participant nucleons, (Np), of about 350 
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A. It was recently observed |T5| that the mean multiphcity of J/ip mesons increases 
proportinally to the mean multiphcity of negatively charged hadrons (more than 90% vr^ 
mesons) when p+p, p+A and A+A collisions are considered (see Fig. 2). It may be 
therefore assumed that a similar dependence on the size of the colliding objects is also 
valid for D-mesons. Following this assumption a mean multiplicity of + D mesons 
can be calculated as: 

{D^ + lf)p,P, = {D' + lf)NN ■ = 2 ■ 10- ■ ^ - 4.5 ■ 10-, (2) 

{n )nn o.l 



where the values of {h )pbPb and {h are taken from Ref. and Ref. respec- 
tively. 

B. It is usually assumed that due to the large mass of the charm quark, production of 
charm can be calculated using perturbative QCD methods |T^. This assumption leads to 



the expectation that the cross section for charm production increases as for all inelastic 
A+A collisions and that the mean multiplicity of open charm increases as (Ap)^/^ for 
central A+A collisions. Thus for central Pb+Pb collisions we get: 

{D'+^)pbPb = {D' + ^)nn ■ ((iVp)/2)^/3 ^ 2 ■ 10-^ (3) 

This estimate agrees with a previously published prediction for charm production based 
on pQCD inspired models ||2^ . 



C. The NA50 Collaboration found recently that a model based on pQCD can not 



describe the dimuon invariant mass spectrum between the and J/ijj peaks in central 
Pb+Pb collisions. The spectrum, however, can be reproduced when the contribution from 
open charm decays is scaled up by a factor of about 3. Based on this analysis one may 
expect that: 

{D' + D')pbPb-Q -10-', (4) 

i.e. it is equal to the multiplicity calculated in B multiplied by a factor of 3. 

D. The production of entropy and strangeness in A+A collisions at the CERN SPS can 
be described by a statistical model which assumes the creation of a Quark Gluon Plasma 
in the early stage of the collision 0]. This fact and the observation that the J/ip yield is 
proportional to the pion yield triggered the hypothesis that also charm production can be 
described by the same statistical approach. It was calculated, within this model p[, that 
the mean number of c and c quarks produced in central Pb+Pb collision at 158 A-GeV is 
about 17. Based on p+p data |T^ we assume here that about one third of them hadronize 
as and mesons which gives: 



+ D")pfep, ^ 6. (5) 



We summarize this section by observing that predictions of the mean multiplicity of 
D'^ + mesons in central Pb+Pb collisions at 158 A-GeV range from 4.5-10"^ to about 
6 showing more than a factor 100 difference between minimum and maximum values. 
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4 D— Meson Measurement in NA49 



We consider here the possibihty of a measurement of and production in A+A 
colhsions at 158 A-GeV using the current set-up of the NA49 experiment at the CERN 
SPS [|[. We start from an introductory discussion on the means to distinguish de- 
cays from combinatorial background and on the calculation of the statistical error on D'^ 
multiplicity. We continue with a brief summary of the relevant properties of the NA49 
experiment and finally we present results of a simulation which yields an estimate of the 
statistical resolution of D^ + D measurement in central Pb+Pb collisions using the NA49 
data. 

Our analysis is based on the study of the two body decay channels of and D 



mesons: 



D° ^K- + 7r+ and D° ^ K+ + vr" (6) 



for which the branching ratio is measured to be (3.85 ± 0.09)% |2^. The mass and the 
proper life time of D° 
cm, respectively p6 . 



proper life time of D° and mesons are m^o = (1864.6 ± 0.5) MeV/c^ and cr = 0.01244 



4.1 Introductory Remarks 

Let us consider central Pb+Pb collision at CERN SPS energy in which among thousands 
of produced particles are also mesons. The mesons decay after a typical flight 
distance of about 0.1 cm (^ 7 ■ cr) from the collision point. There are many decay 
channels possible, but for the reasons discussed at the end of this paper we consider only 
decays of a single type, — * K~ + vr"*", which happen only in 3.85 % of all cases of D° 
decays. 

For simplicity of the initial discussion we assume that an ideal detector is placed around 
the collision point, i.e. for all charged particles electric charge, mass and momentum vector 
at a given reference plane are measured. 

For all {K~, vr"*") pairs originating from decays (signal pairs) the following condi- 
tions have to be approximately (within experimental resolution) fulfilled: 

1. the trajectories of K~ and vr"*" intersect in a point (decay point) which is different 
from the Pb+Pb collision point, 

2. the distribution of the life time in its own c.m.s. frame is an exponential distri- 
bution with a mean value equal to cr = 0.0124 cm, 

3. energy and momentum conservation laws are fulfilled at the decay point when 
decay hypothesis is assumed, 

4. the angular distribution of the decay products in the c.m.s. is isotropic. 

Conditions 1-4 can be fulfilled (within experimental resolution) by a set of {K~ , 7r+) 
pairs which do not originate from the same decay (background pairs) only by chance. 
The probability of this happening decreases with increasing experimental resolution. 
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As a useful example we consider the condition 3 of energy momentum conservation 
at the decay vertex. For practical reasons one often quantifies deviations from energy 
momentum conservation, assuming the decay hypothesis, by calculating the invariant 
mass, M{K~ ,7T~^), for a pair of negatively and positively charged particles and compar- 
ing it to the known mass of the D° meson. Almost all signal pairs are expected to be 
distributed in a narrow interval AM centered around moo ■ The M distribution for back- 
ground pairs is significantly broader and it has no characteristic peak structure at m/jo. 
The smearing of the M distribution for signal pairs results from the finite resolution of 
momentum measurement. With increasing resolution the size of the interval AM in which 
almost all signal pairs are included decreases consequently resulting in lower number of 
accepted background pairs. 

In the high resolution limit almost only signal pairs (and almost no background pairs) 
are selected when taking pairs from AM. In this case the total number of signal pairs, 
Ns, in Nev events can be approximated by the measured number of pairs in the AM 
interval, N^m- An estimate of the mean multiplicity of signal pairs is therefore given by 

(Ni) = ^ (7) 
and its statistical error can be calculated as 

V JV^y V I\ev 

assuming a Poissonian distribution of pair multiplicity in the interval AM. 

In the limit of poor resolution the number of background pairs, Nb, in the interval 
AM is much larger than the number of signal pairs. Therefore in order to estimate the 
number of signal pairs, Ng, an estimate of the number of background pairs, is needed: 

Ni = Nam - N^ (9) 

The estimate of N§ is model dependent and therefore it has a systematic error which 
however will not be discussed here. Various models can be used. As an example let us 
mention a 'mixed event' model of the background in which background pairs are con- 
structed by selecting particles from different events. Independent of the model used it is 
reasonable to assume that the statistical error of N§ can be made much smaller than the 
statistical error of N^m- Thus the statistical error of {Ng) can be calculated as 



V-^^EV y^^EV 

assuming Poissonian distribution of pair multiplicity in the interval AM. 

We observe that in the limit of high resolution the statistical error of the signal is 
ahnost fully defined by the number of signal pairs, whereas in the limit of poor resolution 

it is defined by the number of background pairs. 

One can define the statistical significance of the measurement as the ratio {N§) / a{{Ng) 
which in the case of the poor resolution limit is given by 

(Ni) (Ni) 



■ sJn^. (11) 
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Thus in order to maximize the statistical significance of the result one should select 
the acceptance (in the example above the size of the AM interval) such that the ratio 

{■^s)/\/(^b) reaches a maximum. 



4.2 The NA49 Experiment 

The NA49 experiment at the CERN SPS was designed and constructed to search for 
signals of the Quark Gluon Plasma created at the early stage of nucleus-nucleus colli- 
sions. Basic detectors of the NA49 set-up are four Time Projection Chambers (TPC's), 
which allow for a precise tracking of charged particles. Two medium size TPC's of 3 
m'^ gas volume each are located inside of two magnets with up to 1.5 T field strength 
each. Two large size TPC's of 20 m^ gas volume each are positioned downstream of 
the magnets for high precision energy loss measurement and acceptance coverage in the 
forward direction. Overall the TPC's acceptance coverage amounts to up to 80% of all 
produced charged particles (this number depends on the reaction studied and magnetic 
field - target configuration). 

The TPC's measure up to 234 space points on tracks of up to 13 m length. This 
allows a precise determination of sign of electric charge, momentum vector (cr(p)/p^ ~ 
10~^ (GeV/c)~^) and energy loss (~ 4% relative resolution) for all accepted particles. 
These measurements yield information on particle mass which results in large acceptance 
(limited resolution) particle identification. Four Time-of-Flight walls (~ 60 ps resolution) 
complement the particle identification capabilities of the NA49 detector. 

A typical resolution in reconstruction of the distance between collision point and sec- 
ondary vertex of a two body decay near the target is of the order of 1 cm and is mainly 
due to the long extrapolation length needed from the first TPC detector. 



4.3 Statistical Error on (D^ + D )pupu in NA49 



In order to estimate the statistical resolution of a {D^ + D^) pbpb measurement in NA49 
we performed a simulation of signal and background pairs as expected in this experiment 
for central Pb+Pb collisions at 158 A-GeV. In the simulation the standard geometry and 
magnetic field of NA49 are assumed We define the geometrical accpetance of the 
NA49 TPC's by the requirement that the particle trajectories cross more than 20 TPC 
padrows, at least one of which has to be in a TPC located inside the magnetic field. A 
parametrization of the momentum resolution as measured by NA49 is included in the 
simulation. Concerning particle identification we consider two extreme cases. In the first 
one we assume that no information on particle mass is available, whereas ideal particle 
identification is assumed in the second case. For the background calculation spectra of 
charged hadrons produced in central Pb+Pb collisions at 158 A-GeV as measured by 



NA49 p2| are parametrized in rapidity (y) and transverse momentum (pr)- The signal 



simulation is done assuming a Gaussian {a = 0.6) rapidity distribution and an exponential 



spectrum in transverse mass (T = 300 MeV) for both and mesons. The total 
multiplicity of and mesons of 6 as predicted by the statistical model (see point 
3.2.D) is assumed (if needed). In Fig. 3 we show the y — pr distribution of D° mesons 
for which both decay products are in the geometrical acceptance of the NA49 TPC's. 
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Before further presentation of results of the simulation we note that conditions 1 and 

2 (see Section 4.1), which require an accurate measurement of the decay vertex, can not 
be used in the exisiting NA49 set-up to distinguish signal from background pairs. This 
is due to the poor resolution in the reconstruction of the secondary vertex from two body 
decays of the order of 1 cm, which is much larger than the typical flight path of D° and 
D of about 1 mm. Therefore background rejection has to be done using conditions 3, 
energy-momentum conservation, and 4, isotropy of the decay, only. We quantify these 
conditions by studing M{K^ ^t{^) (or M{K^,t^~)) distribution and cos0 distribution, 
where is the angle between D and tt meson directions calculated in D-meson cm. 
system. 

The results of the simulation shown below are obtained for the D° decay, without 
using information on particle mass. An improvement of the statistical resolution of charm 
measurement due to the addition of D decays and particle mass information is considered 
at the end of this section. 

In Fig. 4 we show an invariant mass distribution M[K~ ,t:^) for background pairs 
obtained using the geometrical acceptance of NA49. It is seen that the meson mass is 
located in the region of the monotonically decreasing tail of the distribution, well beyond 
the position of its maximum. The M{K~ ^ir^) distributions for signal and background 
pairs in the interval of M around itldo are shown in Figs. 5a and 5c, respectively. For 
pairs from the interval AM = 4 MeV/c^ around rriDO we plot also the cosG distributions 
in Figs. 5b and 5d. 

We observe that in the case of measurement by NA49 we are in the poor resolution 
limit i.e. the statistical resolution of signal measurement is determined by the number 
of background pairs and is independent of the signal multiplicity. It is also clear that 
rejection of pairs with high cos© values should result in an improvement of statistical 
significance of the signal measurement. We estimate that the maximum significance can 
be achieved by accepting pairs with cosG < 0.7 and inside an interval of AM = 4 MeV/c^ 
(background cuts). The mean multiplicity of background pairs for this selection is 



where the subscript noid is used to underline that the number is obtained without using 
information on particle identification. This yields an estimate of the statistical resolution 



The resulting correction factors needed to obtain {D^) from the estimated signal mul- 
tiplicity {Ng) in the acceptance are: 

• wsi? = 26.0, for the branching ratio, 

• wqa = 2.4, for the geometrical accepatnce, 

• tysii = 2.3, for the background cuts. 



{Nb)noid ~ 80, 



(12) 



of {Ni) 




(13) 
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Therefore the statistical resolution of (-D°) can be estimated as: 



'{Nb)noid 1300 

(yNOID[{D )) =WbR-WgA-WbC ■ 7^= - rj-j — (14) 

In the case of ideal particle identification background multiplicity can be reduced by 
about a factor of 10, which results in aiD{{D^)) ~ 400/ ^/NEv■ Therefore for 10^ central 
Pb+Pb colhsions at 158 A-GeV a{{D^)) ^ 0.4-1.3 in comparison to ^ 3 expected 
in the case of the statistical model. Statistical significance of measurement of {D^ + D ) 
is approximately by a factor of \/2 better than significance of separate measurements for 
D° or 

In Fig. 6 we plot aNoiD function of the number 

of central Pb+Pb collisions at 158 A-GeV registered by NA49. The different predictions 
concerning {D'^ + D ) discussed in Section 3 are also indicated in Fig. 6 for comparison. 

Finally in Fig. 7 we show M{K~ ,it~^) distributions obtained using the cos0 cut for 
signal and background pairs for 100 central Pb+Pb collisions. In this case the number 
of generated decays was scaled up by a factor of 100 and therefore the statistical 
significance of the observed signal peak corresponds to the significance expected for 10^ 
events. 



4.4 Discussion 

The experiment NA49 registered up to now about 10^ central Pb+Pb collisions at 158 
A-GeV. As follows from the results of simulation presented in Fig. 5 already the analysis 
of this data should yield significant results concerning open charm production in Pb+Pb 
collisions. This encouraging conclusion is reached mainly due to three factors: 

• there are expectations of a significant enhancement of open charm production in 
A+A collisions, 

• a significant fraction of produced Z^-mesons is covered by the large geometrical 
acceptance of NA49, 

• the good momentum resolution of NA49 allows for a significant reduction of back- 
ground even without reconstructing the D-meson decay vertex. 

The results presented in this paper are obtained by a simple analysis of simulated data. 
This scheme was selected in order to underline the main concepts and to build intuition 
concerning basic ingredients of the problem of Z)-meson measurement in A+A collisions. 
The use of more sophisticated statistical methods of data analysis (see Ref. |^), which 



e.g. include explicitly the statistical errors of the measured momentum vector by using a 
kinematical fit, may lead to an improvement of the achieved resolution. 

We studied also the statistical resolution of open charm measurement using three 
and four body decays of neutral and charged D-mesons. We found that due to high 
combinatorial background the statistical resolution is much lower than in the case of the 
two body decay channel. 
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5 Summary and Conclusions 



The main results presented in this paper can be summarized as follows. 

• The measurement of open charm production in A+A collisions gives a unique pos- 
sibility to test predictions of thermodynamical and pQCD-based approaches in the 
region where both are formally applicable. 

• Various approaches used to estimate the D-meson multiplicity in central Pb+Pb 
collisions at 158 A-GeV give predictions which differ by more than a factor 100. 

• Experimental data already registered by NA49 should allow to obtain a significant 
result on open charm production in Pb+Pb collisions. The analysis can yield the 
first observation of open charm signal or it will lead to an estimate of the upper limit 
of open charm multiplicity which should significantly narrow the range of allowed 
models. 

It is obvious that a significant increase of statistics of A+A collisions at maximum 
SPS energy registered by NA49 is required for continuation of open charm program. This 
can be achieved during the already scheduled heavy ion run at the CERN SPS in the year 
2000 and the possible runs beyond. It is also clear that for a precise measurement of open 
charm production an upgrade of the NA49 set-up by a vertex detector, which allows for 
an accurate reconstruction of D-meson decay vertices, is needed. 
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Table 1 The results on mean multiplicity of D + D mesons produced in p+p and 
p+N interactions. The multiplicty for p+N interactions was estimated from the data 
on p+A interactions (used reactions are listed in the third column) assuming that cross 
section for D-meson production is proportional to A. For detailed discussion of the data 
see text. 



Plab [GeV] 


[GeV] 


Reactions 


(DO + ■ 10^ 


Reference 


200 


19.4 


p+Si 


0.47 ± 0.23 


16 




250 


21.7 


p+(Be,Al,Cu,W) 


3.6 ± 0.9 


17 




360 


26.0 


P+P 


6.2 ± 1.9 


18 




400 


27.4 


P+P 


5.6 ± 0.8 


[p^ 




800 


38.8 


P+P 


6.5 ± 2.4 


2G 
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p+(Be,Au) 


5.2 ± 1.0 
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Figure 1: The mean multiplicity of + D mesons in proton-nucleon interactions as a 
function of collision energy in the cm. system. 



16 



xlO" 



> 




20 40 60 80 100 120 

<Np) 



Figure 2: The ratio of the mean muhiphcities of J /if) mesons and negatively charged 
hadrons for inelastic nucleon-nucleon (square) and inelastic 0+Cu, 0+U, S+U and 
Pb+Pb (circles) interactions at 158 A-GeV plotted as a function of the mean number 
of participant nucleons. For clarity the N+N point is shifted from {Np) = 2 to {Np) = 5. 



17 




1 2 3 4 5 6 



y 

Figure 3: The simulated distribution in transverse momentTim and rapidity of D'^ mesons 
for which both decay products {D^ + tt"'") are reconstructed in the NA49 TPC's. 

The distribution is given in arbitrary units. 
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Figure 4: The simulated invariant mass distribution for all combinations of positively 
and negatively charged hadrons accepted by the NA49 TPC's. Masses of positively and 
negatively charged hadrons are assumed to be kaon and pion masses, respectively. The 
calculation is performed for central Pb+Pb coUisions at 158 A-GeV. The spectrum is 
normalized per event and it is in (MeV/c^)~^ units. The shadowed area indicates the 
invariant mass region around mass further used for detailed analysis. 
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Figure 5: The invariant mass and cos© distributions for decays (plots a and b) and 
background (plots c and d) pairs. The cos© distributions arc plotted for pairs from the 
interval AM = 4 MeV/c^ around m^jo. Shadowed areas indicate regions selected for final 
analysis (see text for details). The spectra are normalized per event, the invariant mass 
spectra are in (MeV/c^)~^ units. 
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Figure 6: The statistical resolution of the measurement of mean multiplicity of + D 
mesons in central Pb+Pb collisions at 158 A-GeV as a function of the number of analyzed 
events. The calculation is performed for the current NA49 set-up assuming no information 
on particle mass (upper solid line) and an ideal particle identification (lower solid line). 
The mean multiplicity of + mesons estimated in four different approaches (A— D, 
see Section 3) is indicated by dashed lines. 
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Figure 7: The invariant mass distribution for a sum of signal {D^ decays) and background 
pairs simulated for 100 central Pb+Pb collisions at 158 A-GeV (plots a and b). In the 
simulation the current NA49 set-up was used and the number of signal pairs expected in 
the statistical model (see 3.2.D) was increased by a factor of 100. Thus the statistical 
significance of the result corresponds to the significance expected for 10^ central Pb+Pb 
collisions. The difference of the distribution for a sum of signal and background pairs 
(shown in plot b) and the distribution for background pairs calculated within the 'mixed 
event' model of background is shown in plot c. For all plots only pairs with cos© < 0.7 
are selected. All spectra are normalized per event and they are given in (MeV/c^)~^ units. 
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